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Abstract. Diffuse emission is produced in energetic cosmic ray (CR) interactions, mainly protons and electrons, with the 
interstellar gas and radiation field and contains the information about particle spectra in distant regions of the Galaxy. It may 
also contain information about exotic processes such as dark matter annihilation, black hole evaporation etc. A model of 
the diffuse emission is important for determination of the source positions and spectra. Calculation of the Galactic diffuse 
continuum y-ray emission requires a model for CR propagation as the first step. Such a model is based on theory of particle 
transport in the interstellar medium as well as on many kinds of data provided by different experiments in Astrophysics and 
Particle and Nuclear Physics. Such data include: secondary particle and isotopic production cross sections, total interaction 
nuclear cross sections and lifetimes of radioactive species, gas mass calibrations and gas distribution in the Galaxy (H2, H I, 
H II), interstellar radiation field, CR source distribution and particle spectra at the sources, magnetic field, energy losses, 
y-ray and synchrotron production mechanisms, and many other issues. We are continuously improving the GALPROP model 
and the code to keep up with a flow of new data. Improvement in any field may affect the Galactic diffuse continuum y-ray 
emission model used as a background model by the GLAST LAT instrument. Here we report about the latest improvements 
of the GALPROP and the diffuse emission model. 
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DISCUSSION AND RESULTS 

We give a very brief summary of GALPROP; for details we refer to the relevant papers fU-J^l] and a dedicated 
website. The propagation equation is solved numerically on a spatial grid, either in 2D with cylindrical symmetry in 
the Galaxy or in full 3D. The boundaries of the model in radius and height, and the grid spacing, are user-definable. 
Parameters for all processes in the propagation equation can be controlled on input. The distribution of CR sources can 
be freely chosen, typically to represent supernova remnants. Source spectral shape and isotopic composition (relative 
to protons) are input parameters. Cross-sections are based on extensive compilations and parameterizations [7]. The 
numerical solution is evolved forward in time until a steady-state is reached; a time-dependent solution is also an 
option. Starting with the heaviest primary nucleus considered (e.g., 64 Ni) the propagation solution is used to compute 
the source term for its spallation products, which are then propagated in turn, and so on down to protons, secondary 
electrons and positrons, and antiprotons. In this way secondaries, tertiaries, etc., are included. Primary electrons are 
treated separately. The proton, helium, and electron spectra are normalized to data; all other isotopes are determined 
by the source composition and propagation, y-rays and synchrotron emission are computed using interstellar gas data 
(for pion-decay and bremsstrahlung) and the interstellar radiation field (ISRF) model (for inverse Compton). The 
computing resources required by GALPROP are moderate by current standards. 
Recent extensions to GALPROP include 

• new detailed calculation of the ISRF JHH] 

• proper implementation of the anisotropic inverse Compton scattering using new ISRF (Figure 1, left) 

• interstellar gas distributions based on current HI and CO surveys iflolfTlll 

• new parameterization of the n° production in pp-collisions lfl2tl which includes the diffraction dissociation 

• non-linear MHD wave - particle interactions (wave damping) [6j are included as an option 

• the kinetic energy range is now extended down to ~1 keV 
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FIGURE 1. Left: The ratio of anisotropic IC to isotropic IC for Galactic longitudes 1 = 0° and 180° vs. Galactic latitude. Right: 
y-ray spectrum of inner Galaxy (330° < I < 30°, \b\ < 5°) for an optimized model. Vertical bars: COMPTEL and EGRET data, 
heavy solid line: total calculated flux. This is an update of the spectrum shown in 



• the y-ray calculations extend from keV to tens of TeV (e.g., Figure 1, right), and produce full sky maps as a 
function of energy; the output is in the FITS-format 

• gas mass calibration (Xco-factors) which can vary with position 

• a dark matter package to allow for propagation of the WIMP annihilation products and calculation of the 
corresponding synchrotron and y-ray skymaps 

• GALPROP-DarkSUSY interface (together with T. Baltz) will become publicly available soon 

• a dedicated website has been developed ( |http : //galprop . Stanford. edu| l 

The GALPROP code 01 was created with the following aims: (i) to enable simultaneous predictions of all relevant 
observations including CR nuclei, electrons and positrons, y-rays and synchrotron radiation, (ii) to overcome the 
limitations of analytical and semi-analytical methods, taking advantage of advances in computing power, as CR, y- 
ray and other data become more accurate, (iii) to incorporate current information on Galactic structure and source 
distributions, (iv) to provide a publicly-available code as a basis for further expansion. The first point is the most 
important: all data relating to the same system, the Galaxy, must have an internal consistency. For example, one cannot 
allow a model which fits secondary/primary ratios while not fitting y-rays or not being compatible with the known 
interstellar gas distribution. There are many simultaneous constraints, and to find one model satisfying all of them is a 
challenge, which in fact has not been met up to now. Upcoming missions will benefit: GALPROP has been adopted as 
the standard for diffuse Galactic y-ray emission for NASA's GLAST y-ray observatory, and is also made use of by the 
ACE, AMS, HEAT and Pamela collaborations. 
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